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Abstract. BeppoSAX observed 8 times the 7-ray bright 
blazar PKS 0528+134 in Feb. and Mar. 1997, during 
a multiwavelength campaign involving EGRET, onboard 
CGRO, and ground based telescopes. The source was in 
its faintest X-ray state observed so far, with an unab- 
sorbed [2-10] keV flux of 2.5 x KT 12 erg cm" 2 s -1 . The 
spectrum can be fitted with a power law, with energy in- 
dex a x = 0.49 ± 0.07 between 0.1 and 10 keV. The high 
energy instrument PDS detected a flux in the 15-100 keV 
band which appears disconnected from the lower X-ray 
energies, most likely due to a contaminating source. We 
consider our findings in the context of the overall spec- 
tral energy distribution and discuss their implications for 
synchrotron and inverse Compton models. 
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1. Introduction 

Blazars received an unexpected boost of interest since the 
discovery, driven by EGRET (and confirmed by OSSE and 
COMPTEL) that they emit the bulk of their radiative 
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power at 7-ray energies. We are now finding spectral dif- 
ferences among blazars which have already been detected 
by EGRET. BL Lac objects seem to be characterized by 
a 7-ray luminosity, L 7 , which is comparable to the lumi- 
nosity emitted in the rest of the spectrum, while L 7 tends 
to dominate in blazars with visible broad emission lines. 
At the same time, the 7-ray dominance correlates with 
the overall spectral energy distribution (SED) (Comastri 
et al. 1997, Fossati et al. 1998). All blazar SED are in 
fact characterized by two broad emission peaks, believed 
to be produced by the synchrotron and the inverse Comp- 
ton processes. The location of these peaks (their frequency 
and relative flux) is a strong diagnostic tool to discrimi- 
nate among theoretical models and to find the intrinsic 
physical parameters of the emitting region (see e.g. Ghis- 
ellini, Maraschi & Dondi 1996; Sikora et al. 1997; Dermer, 
Sturner & Schlikeiser 1997). Examining the SED of 7— ray 
detected blazars, one finds that both peaks are located at 
lower frequencies in more powerful, and more 7-ray dom- 
inated, sources (Fossati et al. 1998; Ghisellini et al. 1998). 

PKS 0528+134, a very powerful blazar, is no exception: 
its SED shows that the high energy peak is located in the 
10-100 MeV range, and the low energy peak between the 
far IR and the optical bands. For this source, the X-ray 
band falls near the minimum between the two peaks. Pre- 
vious X-ray observations by the ASCA satellite showed a 
flat (a < 1, with F v oc v~ a ) spectrum, possibly connect- 
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ing with the flux in the MeV-GeV band (Sambruna et al. 
1997). The uncertainty comes from the paucity of simulta- 
neous observations in both bands, and from the relatively 
narrow energy windows of the used X-ray detectors. In- 
deed, if one could join the X-ray and the 7-ray emission 
with a smooth curve and, more importantly, simultaneous 
variability were detected, then it would be possible to con- 
fidently argue that the emission in both bands is produced 
by the same mechanism and in the same region. 

Furthermore, the fact that the X-ray emission of this 
source is in a minimum of the SED is of particular interest. 
In fact, the possible reprocessing of high energy photons, 
via the creation of electron-positron pairs, leads inevitably 
to the production of radiation, mainly in the X-ray band. 
It is therefore crucial to determine how deep this minimum 
is: a deep minimum would be a signature that this process 
does not occur (Ghisellini & Madau 1996). 

For these reasons we have undertaken a program 
to observe PKS 0528+134 with the X-ray broad band 
BeppoSAX satellite simultaneously with CGRO (in par- 
ticular with the EGRET instrument) and other, ground- 
based, telescopes. In this paper we mainly present the 
BeppoSAX observations, compare them with previous X— 
ray observations, and construct a simultaneous SED with 
the optical and EGRET data, briefly discussing our find- 
ings. 

2. PKS 0528+134 

PKS 0528+134 (z = 2.07) is one of the most distant 
quasar detected by EGRET in the 7-ray band. Located 
in the galactic anticenter region, it is heavily absorbed, 
although the estimates of Ay are very uncertain (ranging 
from 2.3 to 5). Thus PKS 0528+134 is faint in the optical, 
with a typical average magnitude of 19.5 in the V band. It 
is a strong and flat radio source, with detected superlumi- 
nal motion (/3 app ~4-6 for H = 100 km s" 1 Mpc" 1 , Pohl 
et al. 1996; Britzen, Witzel & Krichbaum 1996). Close 
to Geminga and the Crab, it was frequently observed by 
EGRET, and seen flaring in March 1993, when the 7-ray 
flux was about 6 times brighter than average (Hunter et al. 
1993; Mukherjee et al. 1996). Flux changes in the 7-band 
are accompanied by spectral variations in the sense that 
the spectrum is harder when brighter (Sambruna et al. 
1997; Mukherjee et al. 1997a, 1997b). It is also one of the 
few extragalactic sources detected by COMPTEL (Coil- 
mar et al. 1997; Bottcher & Collmar 1998). The combined 
COMPTEL and EGRET observations show that the peak 
of the high energy emission occurs in the MeV band, at 
least in the flaring state. 

3. BeppoSAX. observations 

Between Feb. 21 and March 11, 1997, BeppoSAX observed 
the source 8 times with the narrow field instruments LECS 
(Low Energy Concentrator Spectrometer; energy range 



0.1-10 keV), MECS (Medium Energy Concentrator Spec- 
trometer; energy range 1.3-10 keV) and PDS (Phoswich 
Detector System; energy range 15-200 keV) (Boella et al. 
1997). 

In Table 1 we report the journal of observations, with 
the corresponding exposure times and derived net count 
rates. 



0.05 



s" 0.045 

60 

S 0.04 

O 

0.035 



17.8 

18 

-" 18.2 

18.4 

18.6 

8 



S ft P 



■S 4 



MECS 
Net cot 



unt rate 1 ,, 



I I I I I I I I I I I I I I I I I I I I 
Optical T 1 

M } S 1 ; 

I I I I I I I I I I I I I I M I I I I 



, Radio 



8.3 GHz 



2.25 GHz 

I 1 1 1 1 I 



eeee eeeeg 



_L 



_L 



eeee e_J 



10 20 30 40 

Days of Feb. 1997 



50 



Fig. 1. The top panel shows the MECS light curve, corre- 
sponding to the 8 observations. In the middle panel we report 
the optical light curve in the same period, in the Ic filter. The 
bottom panel shows the light curve at 8.3 and 2.25 GHz, by 
the Green Bank monitoring campaign. 



4. Results 

4-.1. Light curves 

We have extracted light curves from the data of the three 
MECS units (merged together) using XSELECT and con- 
sidering a circle of radius 4' centered around the target. 
In Fig. 1 we show the total MECS count rates obtained 
with the XIMAGE program corresponding to the 8 ob- 
servations. This light curve is consistent with a constant. 
We can set an upper limit on the amplitude of possible 
variations of ~40%. No significant variations are present 
within the single observations either. 

Simultaneously with the BeppoSAX observations, the 
source was observed in the optical, with the REOSC 1.05 
meter telescope of Pino Torinese (see Villata et al. 1997 
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TABLE 1: BeppoSAX observation log 



Date 


Start 1 


End 1 


LECS^ 


net cts/s 


MECS^ 


net cts/s 


PDS X 


net cts/s a 


Feb 21 


02:17:38 


09:08:06 


6322 


0.015±0.007 


14487 


0.042±0.002 


6919 


0.22±0.08 


Feb 22 


23:29:16 


06:48:06 


5046 


0.014±0.007 


13365 


0.045±0.003 


5925 


0.24±0.08 


Feb 27 


14:38:39 


19:18:11 


3945 


0.011+0.010 


6917 


0.041±0.004 


1622 


0.28±0.17 


Mar 1 


00:52:14 


07:43:06 


4261 


0.012±0.010 


13547 


0.043±0.003 


6344 


0.29±0.08 


Mar 3 


01:38:02 


08:28:06 


5053 


0.015±0.007 


14241 


0.039±0.002 


6783 


0.16±0.08 


Mar 4-5 


21:53:08 


04:28:05 


2793 


0.015±0.010 


11254 


0.043±0.003 


4914 


0.11±0.09 


Mar 6-7 


21:03:14 


03:53:05 


3083 


0.012±0.010 


12762 


0.041±0.003 


5587 


0.26±0.09 


Mar 11 


00:17:06 


06:06:05 


2497 


0.012±0.011 


11490 


0.037±0.003 


5374 


0.29±0.09 


Total 






32998 


0.014±0.001 


98063 


0.041±0.001 


43468 


0.23±0.03 


1 UT Time; 2 Exposure time 


in seconds 


a values in the entire PDS band 







for a description of the instrument and the data analysis 
procedure), and in the radio band with the Green Bank 
Interferometer (GBI). As shown by the middle and lower 
panels of Fig. 1 , the source was varying significantly in the 
optical (with a probability > 99.9%) and at 8.3 GHz. 

We also observed the source in the near infrared with 
the 1.9 m telescope at Sutherland using the Mklll pho- 
tometer, obtaining K = 14.76 ± 0.23 on Mar 2, 1998 
(JD=2450510), using a 12" diameter aperture. 

4-2. Spectral fitting 

4.2.1. LECS and MECS 

Spectra of the LECS and MECS detectors have been ex- 
tracted using XSELECT and a circle of radius 4' around 
the target for both the instruments (given the weak flux 
of the source and the heavy absorption under 0.7 keV, we 
have not used the 8.5' radius recommended for the LECS 
for brighter sources, to minimize the effects of the back- 
ground, especially at low energies). As suggested by the 
Science Data Center of BeppoSAX (SDC), after checking 
that the background was constant during the observations, 
we have subtracted the background using the blank sky 
fields selecting the same detector regions where the target 
is located. We used the response matrices released by the 
SDC in September 1997. The LECS and MECS spectra 
have been jointly fit after allowing for a constant rescaling 
factor of 0.9 for the LECS data (this factor accounts for 
uncertainties in the inter-calibration of the instruments; 
Fiore, priv. comm.). 

Fitting the 8 observations separately with a single 
power law and free hydrogen column density, Ah, gives 
consistent values, within the (rather large) uncertainties. 
We therefore decided to combine the 8 spectra, in order 
to increase the signal to noise ratio. A single power law fit 
to the LECS+MECS dataset (see Fig. 2) yields an energy 
spectral index a x = 0.48±0.11 with A H = (5.0±1.9)xl0 21 
cm~ 2 (unless otherwise indicated, the errors from the fits 
are at 90% confidence level, for two parameters of inter- 
est). No spectral feature is required by the fit and in par- 
ticular an upper limit of EW < 200 eV can be set to the 



equivalent width (in the source frame) of the 6.4 keV flu- 
orescence K a iron line (see §5.1 for more discussion). The 
results of the spectral fitting are reported in Table 2. 



PKS 0528+134 
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Fig. 2. Fit to the LECS and MECS BeppoSAX spectrum. 
PDS data are not used for the fit, and appear to lie above the 
extrapolation of the power law from lower energies. 



As mentioned in the Introduction, the value of the 
Ah column density is uncertain, yet its determination is 
important, since the knowledge of the corresponding op- 
tical extinction allows to determine the location of the 
synchrotron peak. BeppoSAX data, alone, cannot well 
constrain Ah, due to the faint level of the source. PKS 
0528+134 was instead much brighter during the ASCA 
observations of March 1995, allowing a better determina- 
tion of the A H column [A H = 5.3(±0.3) x 10 21 cm" 2 ], as 
described in §5.1. 



Ghisellini et al., BeppoSAX observations of PKS 0528+134 



TABLE 2: Fits to BeppoSAX Data: single power law 



7V H F { 

10 21 cm -2 erg/cm 2 /s 



[2-10 keV] Xr 



d.o.f. 



■Fl koV 

fiJy 



Notes 



0.48±0.11 5.0±1.9 

0.49±0.07 5.3 

0.37±0.11 3.9±1.6 

0.42±0.06 5.3 



2.53e-12 


0.91 


54 


2.54e-12 


0.91 


55 


2.59e-12 


1.71 


58 


2.60e-12 


1.74 


59 



0.29±0.05 LECS+MECS, free N H 

0.30±0.03 LECS+MECS, fixed N H 

0.24+0.04 LECS+MECS+PDS, free Nu 

0.27+0.03 LECS+MECS+PDS, fixed N K 



TABLE 3: Single power law fits to ASCA data with fixed JV H 



Q 


F[2-10 kcV] 

erg cm' 2 s _1 


xi 


d.o.f. 


-Fl kcV 

fiJy 


Date 


0.66+0.06 


3.77e-12 


0.89 


93 


0.57+0.05 


28 Aug 1994 


0.71+0.04 


9.55e-12 


1.21 


100 


1.57+0.07 


7 Mar 1995 


0.71+0.04 


10.62e-12 


1.18 


69 


1.76+0.07 


14 Mar 1995 


0.72+0.03 


13.80e-12 


0.91 


140 


2.32+0.07 


19 Mar 1995 



Fits to the summed ASCA data of March 95 



<U 



"2 



-E^break 

keV 



10 21 c 



F, 



[2-10 kcV] Xr 

erg/cm 2 /s 



d.o.f. Notes 



0.72+0.04 - 
0.00+0.16 0.65+0.03 



1.47 + 0.10 



5.3+0.3 1.15e-ll 1.13 207 free JV H , single power law 

3.9 fix 1.15e-ll 1.08 206 fixed N H , broken power law 



Since the BeppoSAX measurement of the Nu column 
is consistent with the value found by ASCA, we fixed the 
TVh at the latter value. In this case the LECS+MECS data 
can be modeled by a single power law with an energy index 
a x = 0.49 ± 0.07 (90% confidence level for one parameter 
of interest). 

4.2.2. PDS 

Source visibility windows were selected following the cri- 
teria of no earth occultation and high voltages stability 
during the exposures. In addition, the observations clos- 
est to the South Atlantic anomaly were discarded from the 
analysis. From collimators positions the ON and OFF time 
windows were also created and merged with the source vis- 
ibility window to create the final time windows on which 
the source+background and background spectra were ac- 
cumulated for each of the four PDS units, using the XAS 
software package. A filtering for the temperature and en- 
ergy dependence of the pulse rise time was used. 

The source was significantly detected up to ~90 keV. 
For each pointing, the grouped spectra from the four units 
were coadded. The net count rates are reported in Table 
1. 

To increase the S/N, we have averaged the eight PDS 
spectra and binned the resulting spectrum in 4 energy 
intervals. For the fitting, we used the response matrix 
available at SDC, keeping the rescaling factor with re- 
spect to the MECS constant and equal to 0.85 (this factor 
allows for uncertainties in the inter-calibration of the in- 
struments). 



As shown in Fig. 2, the PDS data points lie above the 
model fit for the LECS and MECS datasets. In particular, 
including the PDS data and using a single power law fit 
with iVn fixed to the value found by ASCA, we obtain 
a x = 0.42 ± 0.06, with the PDS points still lying above 
the fit (see Table 2 for details). 

The same three datasets (LECS+MECS+PDS) were 
also fitted with a broken power law model with a\ = 0.48+ 
0.06 and a 2 = —0.15 ± 0.14 and break energy Eb re ak ~ 
9i| keV, leading to a xl = 1-34 for 57 d.o.f. This model 
significantly improves the fit, according to the F-test (> 
99.5%), but the x 2 ~value yields a probability less than 5% 
that the model is correct. A power law fit to the PDS data 
alone gives a spectral index a x — 0.82 ± 0.35. 

As can be seen, the high energy flux detected by 
the PDS shows a clear discontinuity with respect to the 
LECS/MECS data. A possible explanation can be that 
the PDS data are contaminated by the presence of an- 
other source in its field of view. The PDS field of view 
is ~1.3 degree, with no imaging capability. The response 
matrix of this instrument is triangular, with a flat top of 
3' and a reduction of a factor 2 in sensitivity at 38' from 
the center (Frontera et al. 1997). The MECS have a field 
of view of radius 30', and there are no sources other than 
PKS 0528+134 detected at more than 3cr. 

We checked the ROSAT PSPC WGA catalog (White, 
Giommi & Angelini, 1994) for sources at an angular dis- 
tance < 1 degree from PKS 0528+134. There are two 
unidentified sources (1WGA J0532. 9+134, at a distance of 
31 arcmin from PKS 0528+134, and 1WGA J0533+135, 
at 41 arcmin) with a count rate in the ROSAT band com- 
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parable to that of PKS 0528+134. No sources were instead 
found in the EXOSAT and ASCA public archives. Given 
the low galactic latitude of PKS 0528+134 and the large 
offset of the PDS data points with respect to the extrapo- 
lation of the MECS spectrum, we think that the PDS data 
might be contaminated by other sources. For this reason, 
in the discussion of the SED of PKS 0528+134, we will 
not consider the PDS spectrum. 

5. Comparison with previous observations 

5.1. ASCA observations 

As already mentioned, a long standing problem with this 
source is the determination of its optical extinction and 
the column Ah along the line of sight. The value that 
we found with the BeppoSAX. observations is in agree- 
ment with the sum of the estimated absorption caused 
by the column of neutral galactic hydrogen [Nu — (2.6 ± 
0.1) x 10 21 cm -2 ], plus that due to the outer edge of the 
molecular cloud Barnard 30 in the A Orion rings of clouds, 
[N H ~ 1.3 x 10 21 cm -2 ] (Sambruna et al. 1997). 



cm -2 ], without considering the extra absorption due to 
the Barnard cloud. 

We therefore decided to re-analyze the full ASCA data 
set. We first fitted the 3 observations separately, find- 
ing that, despite a flux change of a factor 1.5, the three 
best fit values of the spectral index were nearly equal: 
with free Ah the values were within the 68% confidence 
uncertainties. We then fitted the sum of the 3 ASCA 
observations and the increased signal-to-noise ratio al- 
lowed to obtain the best measure of the column so far: 
A H = (5.3 ± 0.3) x 10 21 cm -2 (see Table 3 for the fit 
parameters). 

In Fig. 3 we show the confidence contours in the pho- 
ton index-column plane, for the sum of the August 1995 
ASCA observations and the sum of the 8 BeppoSAX ob- 
servations. 

We also searched for the possible presence of an Fe 
K a line at 6.4 keV, finding an upper limit to its equiva- 
lent width of 67 cV in the summed 1995 spectrum (calcu- 
lated in the rest frame of the source). This limit is only 
marginally consistent with the value of EW = 119 + 58 
eV reported by Reeves et al. (1997) analyzing the 1994 
ASCA spectrum, when the source was fainter. 
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Fig. 3. Confidence contours (68%, 90% and 99%) for fits with 
a single power law plus absorption for the sum of the three 
August 1995 ASCA observations and the sum of the eight 
BeppoSAX observations. 



In March 1995, ASCA observed PKS 0528+134 three 
times, within 2 weeks, during an active state of the source. 
The last of these observations was analyzed by Sam- 
bruna et al. (1997), finding a best fit value of Ah = 
(5.57 ± 0.4) x 10 21 cm -2 , in good agreement with the 
value found in our BeppoSAX observations. The entire set 
of the three observations was instead analyzed by Kubo 
et al. (1998), but assuming a column density fixed at the 
neutral galactic hydrogen value [Ah = (2.6 ± 0.1) x 10 21 



5.2. Historical X-ray and j-ray light curves 

In order to understand the relationship between X-ray 
and 7— ray emission, we compared both the light curves 
in the last 7 years and the spectral indices behaviour in 
these two bands. 

In Fig. 4 (upper panel) we show the light curves of all 
the observations in the X-ray band after 1991. For con- 
sistency, we have re-analyzed all these data (taken from 
the archives) with the same Ah = 5.3 x 10 21 cm -2 . In 
1995 the source varied by 50% in 2 weeks. Between the 
1995 ASCA and our 1997 BeppoSAX observations, when 
the source was at its faintest historical level, the flux de- 
creased by a factor ~7. The lower panel shows the light 
curve of the 7-ray flux, as observed by EGRET (Mukher- 
jee et al. 1996, 1997a, 1997b). In 7 rays the source was 
observed more often, with a variability of a factor 13 in 2 
months, and a factor ^2 in 2 days during the 1993 flare 
(see Mukherjee et al. 1996). 

The lack of simultaneity of previous observations and 
the paucity of X-ray data do not allow to establish any (or 
absence of) correlation. It should be however noticed that 
for the three epochs of simultaneous observations (Aug. 
1994, Mar. 1995 and Feb.-Mar. 1997), the X-ray and 7- 
ray flux levels follow the same trend. 

5.3. Flux-spectral index correlations? 

In Fig. 5 (upper panel) we report the X-ray spectral in- 
dices determined by BeppoSAX and ASCA data vs the 
1 keV flux. While all ASCA points are consistent with 
the same spectral index, the spectral shape determined 



Ghisellini et al., BeppoSAX observations of PKS 0528+134 




Time 



Fig. 4. Historical light curves in the X-ray and 7-ray bands, 
after 1991. All X-ray data have been reanalysed assuming a 
fixed JV H = 5.3 x 10 21 cm -2 . 
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by BeppoSAX observations is significantly flatter (see also 
Fig. 3), and corresponds to the fainter flux. 

In the lower panel of Fig. 5 we show the 7-ray spectral 
index as a function of the 7-ray flux (data from Mukherjee 
et al. 1996, 1997a, 1997b, 1999). The large error bars of 
a 7 , especially at faint fluxes, do not allow any conclusion 
about the presence of a trend: although the linear cor- 
relation coefficient is ~-0.6 (random probability =0.02), a 
X 2 -test on a fit with a constant gives 15% probability that 
the two quantities are not correlated. 

Although we cannot draw any firm conclusion about 
the flux-spectral index correlations, note that during the 
BeppoSAX observations the source, in its faintest state, 
had the flattest spectral index. This is an unusual behav- 
ior for blazars (even if some other example exist, see Ul- 
rich, Maraschi & Urry 1997), and likely to yield important 
information and/or constraints on the emission models, as 
discussed below. 

6. The SED 

In Fig. 6 we show the overall spectral energy distri- 
butions of PKS 0528+134 at different epochs, as indi- 
cated by the labels. Data have been dereddened assuming 
Nn = 5.3 x 10 21 cm~ 2 , corresponding to Ay — 3. We have 
reported the minimum and maximum values of the radio 
and optical data obtained during our campaign. 

The not simultaneous data shown in Fig. 6 have been 
taken from Wall & Peacock 1985; Bloom et al. 1994; Edel- 
son 1987; Reuter et al. 1997; Rieke et al. 1982; Condon et 



Fig. 5. Spectral indices vs flux in the [2-10] keV and [0.05-1] 
GeV bands. All X-ray data have been reanalysed assuming a 
fixed Nn = 5.3 x 10 21 cm -2 . 7-ray data from Mukherjee et al. 
1996, 1997a, 1997b, 1999. 



al. 1977; Collmar et al. 1997; McNaron-Brown et al., 1995 
and Mukherjee et al. 1996. We also show some simultane- 
ous sets of data: the 1991 Comptel and EGRET spectrum 
(Collmar et al. 1997); the 1994 and 1995 X-ray (ASCA) 
and EGRET data (we also have an optical point during 
the 1994 campaign; data for these two campaigns are pre- 
sented in Sambruna et al. 1997 and references therein); our 
1997 IR, optical, X-ray data together with the EGRET 
point from Mukherjee et al. 1997b and the radio data from 
the GBI archive (see Fig. 1). 

As all other 7-ray bright blazars, also the SED of 
PKS 0528+134 is characterized by two peaks, one between 
the far IR and optical spectral bands, and the other at 
MeV energies. As discussed in Sambruna et al. (1997), 
the dereddened (using Ay = 3) optical spectrum is in- 
verted (a — —0.18 ± 0.08), and thus indicative of the 
presence of a 'blue bump' component. This is not easy to 
reconcile with the very large and rapid optical variabil- 
ity. We also note that due to the presence of molecular 
clouds along the line of sight the gas to dust ratio may be 
anomalous and the value of Ay may be larger than the 
one adopted above. Further spectral observations in the 
IR optical range could help clarifying this issue. 

As can be seen, the MeV-GeV emission dominates 
the bolometric output by a large amount, reaching an 
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Fig. 6. Overall spectral energy distribution of PKS 0528+134. Simultaneous observations are indicated by the different marks, 
as labeled. Other non-simultaneous data are taken from the literature (see text). The solid line is a fit to the 1997 spectrum 
made with a homogeneous, one-zone model, corresponding to a source with a size R — 5 x 10 16 cm, a magnetic field B — 8.6 
Gauss, in which beaming corresponds to a Doppler factor S = 17. Relativistic electrons are injected continuously throughout the 
source, at a rate corresponding to an intrinsic luminosity L' = 3.7 x 10 44 erg s , with a power law spectrum oc 7~ 3,8 between 
71 = 200 and 72 = 3000. fn this case the EC component dominates the entire high energy band. The fit to the 1995 spectrum is 
obtained adopting the same R, 8 and B, but with a larger injected power of V = 1.8 x 10 45 erg s _1 and with a flatter injected 
electron spectrum (oc 7~ 2 ' 7 ) between 71 = 200 and 72 = 4000. fn this case both the SSC and the EC components contribute to 
the low energy X-ray spectrum, while only the EC component contributes above a few tens of keV. Note that the shown PDS 
points lie above the extrapolation from the spectrum in the LECS and MECS band. Given the possibility of contaminating 
sources (see text), we do not consider these points when applying the model. 



(isotropic) luminosity in excess 
Hq = 50 km s _1 Mpc -1 and qo = 



7. Discussion 



of 10 49 
0.5). 



erg s _1 (with & Celotti 1992; Bloom & Marscher 1993; Sikora, Begel- 
man, & Rees 1994, Dermer, Schlickeiser & Mastichiadis 
1992; Blandford & Levinson 1995; Ghisellini & Madau 
1996). 



Among the models proposed to account for the blazar 7- 
ray emission, the most accredited are those based on the 
Inverse Compton process, between the relativistic elec- 
trons responsible for the synchrotron emission and seed 
photons which can cither be the synchrotron photons 
themselves [synchrotron-self Compton (SSC)] or photons 
produced externally to the jet (EC) (Maraschi, Ghisellini 



In all cases synchrotron radiation is responsible for the 
low energy emission (up to the optical-UV band), and the 
inverse Compton process gives rise to the X- to 7-ray spec- 
trum. Seed photons of different origin may well coexist 
and their relative contribution in the same object may be 
different in different states and in different energy bands. 
The presence of luminous emission lines and possibly of a 
blue bump indicates, in the case of PKS 0528+134, that 
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photons outside the jet should play an important role in 
the IC process (EC scenario). For a detailed discussion of 
the fitting of the 1994 overall spectrum with SSC and EC 
models see Sambruna et al. (1997). 

As shown in Fig. 6, the source shows correlated X-ray 
and 7-ray flux variability. In particular, during the 1997 
campaign, the source was at its faintest level both in X- 
and 7-rays. In this occasion, the X-ray spectrum was sig- 
nificantly flatter than in higher states. This is in contrast 
with the "flatter when brighter" behavior of blazars (Ul- 
rich, Maraschi & Urry, 1997), and therefore likely to give 
important constraints on the emission models. 

Within the framework of the SSC and EC model we 
can envisage two scenarios to explain the observed behav- 
ior: 



1) The X-rays are produced by both the self-Compton 
and EC processes. Therefore there can be two typical fre- 
quencies of the seed photons: one corresponding to the 
peak of the synchrotron emission (in the far IR), and one 
corresponding to the external radiation. In the case that 
the emission lines and blue bump photons form the bulk 
of the external radiation, their typical frequency, as ob- 
served in the comoving frame, is in the far UV. If the ra- 
diation energy densities of these two components are com- 
parable (within a factor 10), then the self-Compton emis- 
sion would dominate at lower X-ray energies, while the 
EC spectrum would be entirely responsible for the emis- 
sion in the 7-ray band. This is because the self-Compton 
spectrum is somewhat steeper than the EC one produced 
by the same electrons, and because the maximum self 
Compton frequency is lower. Let us suppose now that the 
number of emitting electrons increases. In this case the 
synchrotron and the EC fluxes vary linearly, while the 
self-Compton one varies quadratically, making the SSC 
component to dominate the flux over a larger X-ray en- 
ergy range. 

We have reproduced the 1997 "low state" and the 1995 
"high state" spectra along these lines, with a homoge- 
neous one-zone model, taking into account both the SSC 
and the EC contributions to the high energy spectra. The 
model finds the equilibrium (steady state) electron distri- 
bution by balancing the injection of relativistic particles 
and their radiative cooling. Photon-photon collisions and 
electron-positron pair production, Klein Nishina effects 
on the scattering cross section and Coulomb collisions are 
taken into account (see Sambruna et al. 1997 for further 
details). The input parameters are given in the caption of 
Fig. 6. We "fitted" both spectra leaving the dimension and 
the beaming factor of the source unchanged, as well the 
low energy cut-off of the injected electron distribution. We 
have instead increased by a factor 3 the injected power in 
the high state, also characterized by a slightly flatter in- 
jected electron distribution. The resulting spectra, shown 



in Fig. 6, reproduce quite well the different X and 7-ray 
slopes. 

2) Another process that can in principle account for a 
flatter X-ray spectral index when the source is fainter, is 
e pair production. When the source is in a high (and 
hard) 7-ray state, a small fraction of the high energy 
power could be converted in electron-positron pairs via 
7 — 7 interaction with X-ray photons. This process pro- 
duces a steepening of the electron spectrum towards lower 
energies, yielding a brighter and steeper IC spectrum in 
the X-ray band. The main difficulty with this scenario is 
the fine tuning required in the amount of the high energy 
power that is absorbed and reprocessed, which must be of 
the order of a few per cent. If it is less, then pairs are not 
produced in a sufficient amount to contribute anywhere 
in the emitted spectrum. If it is higher, then they could 
increase the number of target photons for 7-7 collisions, 
inducing a catastrophic cascade, and an overproduction of 
X-ray with respect to what we observe (Svensson 1987). 

For these reasons, we conclude that the required mod- 
erate pair reprocessing is unlikely to occur in this source. 

8. Summary and Conclusions 

From the analysis of the BeppoSAX and ASCA data, and 
from the construction of the simultaneous SED of the 
source at different epochs we find: 

- BeppoSAX observed PKS 0528+134 in its faintest X- 
ray state, characterized by the flattest X-ray spec- 
trum, of energy index a x = 0.49 ± 0.07. 

— The three ASCA observations of March 1995, despite 
the different fluxes, are characterized by the same 
spectral index a x = 0.71 ± 0.04. Fitting the sum of 
these data allows to determine with improved precision 
the absorbing hydrogen column density of the source, 
N-R = (5.3 ± 0.3) x 10 21 cm -2 , which corresponds to 
an optical absorption of Ay — 3 ± 0.17. 

— The overall simultaneous SED follows the general 
trend of blazars, as found by Fossati et al. (1998). In 
fact the SED is characterized by two peaks, in the mid 
or far infrared and in the MeV band, respectively, as 
for the most luminous blazars. Furthermore the 7-ray 
luminosity dominates the overall power output during 
active states, again in agreement with the general be- 
havior of other powerful blazars (see e.g. von Montigny 
et al., 1995). 

- We have fitted the 1995 SED (high state) and the 1997 
SED (low state) with a homogeneous model, in which 
the emission from the far IR to the UV is produced 
by the synchrotron process, and the high energy flux 
is due to the sum of the synchrotron self-Compton 
and the external Compton contributions. The former 
becomes more important (in the X-ray band) for in- 
creasing source power. This accounts for the unusual 
observed X-ray behavior (flatter state when the source 
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was fainter). Variations of the 7-ray flux and spectral 
shape might be explained by a flattening of the high 
energy part of the electron distribution. 
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